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OF  SHOftT  PULSE 
IN  PLASMAS 


m 


of  Intense  radiation  boost  with  plaaaas  ploys  an 
rolo  In  eooh  sross  ss  lssor  Or  Ivon  fusion,  ionosphorio  prooossos  and 
pill Sir  rsitsTtnn  ftaeently  ths  possibility  of  eaploying  intonss  lssor  beams 

«•  —osier  ale  olootrons  to  ultra-high  snort  its  has  also  stlnulatod  intoroat  In 
radiation  self -focusing  in  plaenas.1*1*  Ono  such  lasor  dr i von  aoooloration 
is  ths  lssor  boot  wave  aooolorator.  This  la  a  collective  aoooloration 
wttloh  utillsos  a  largo  amplitude  plaaaa  wave  with  phase  voloolty 
slightly  loss  than  thO  voloolty  of  light  to  aooolorato  charged  part  idea.  The 
largo  aoplitudo  plaaaa  wavs  is  exoited  by  the  nonlinear  coupling  of  two 
intense  lssor  bosas  propagating  through  the  plaaaa.  In  this  process  the  two 
lasor  boons  with  frequencies  «*1 ,  *2  and  corresponding  wave  numbers  k1 ,  k2 
couple  through  the  plaaaa  to  produce  a  ponderoaotlve  wavo  with  frequency 


t*1  *  *2  11114  oavo  number  k1  -  kg.  If  ^ 


the  plaaaa  wave  will 


Initially  grow  linearly  in  tli 


If  the  laser  frequencies  are  much  greater 


than  the  ambient  plaaaa  frequency,  then  the  phase  velocity  of  the 
ponderoaotlve  wave  la  nearly  equal  to  the  group  velocity  of  the  laser  wave. 
Eleotrona  which  are  either  injected  into  the  plaaaa  or  part  of  the  thermal 
tall  of  the  plaaaa  distribution  can  be  accelerated  by  the  large  gradients 
associated  with  the  plasma  wave. 

In  this  paper  the  self ‘focusing  properties  of  a  plasma  on  an  intense 
electromagnetic  beam  are  studied. *5*22  propagation  of  an  Intense 

radiation  beam  in  a  plasma  can  result  in  self-focusing  by  creating  a  density 
depression  in  the  plasma  as  well  as  by  Increasing  the  electron  mass  by 
relativistic  effects.  The  density  depression  is  due  to  transverse 
ponderomotive  forces  which  tend  to  expel  plasma  from  high  field  regions.  The 
radial  ponderomotive  force  on  the  electrons  is  larger  than  on  the  ions  by  the 
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">  w* qtyMtMNi  farmed  out  of  tt*  radiation  boom  region  oot  up  on 

'■#  • 

rWlrti  faro#  onion,  oo  o  slower  tint  ooolo,  oouooo  tho  ions 
Tho  anon—  ooporotloo  onion  ooto  up  tho  eleetrootatlo 
mtmm  »rt>Xi  mr  o  M^o  looctn  ohloh  lo  ssauusd  saall  ooapared  t< 
Mo  rmUU  W  noon  mdluo.  Thlo  density  doprooolon  orootoo  o  loool  increase 
la  tho  effootivw  index  of  rofr notion  and  ooto  oo  on  optiool  guide  for  tho 


PMttoUoi  noon,  la  addition  to  thlo  self -focusing  aoehonion  o  furthor 
raduatlon  la  tho  ploona  frequency  ooeuro  in  regiono  of  high  fiold  lntonolty 
duo  to  tho  rolotlviatio  aoao  lnoroooo  of  tho  olootrono  in  tho  proaonco  of  tho 
radiation  boon.  Tho  self -focusing  duo  to  tho  donolty  doprooolon  oooura  on  o 
loagor  tlao  ooolo  than  dooa  tho  relatlviotlo  aooo  lnoroooo  self-foeusing 
of  foot.  When  tho  pondoroootlve  proaouro  lo  auoh  loao  than  tho  plaoao  klnotlo 


prooouro  tho  tlao  ooolo  for  tho  donolty  doprooolon  to  occur  lo  tg  •  R/Cg, 
shore  R  lo  tho  radiation  boon  radius  and  Cs  la  tho  Ion  ooouotle  spood.  Slnoo 


R  lo 


auoh  grootor  than  tho  Dobyo  length,  tho  oharaotorlatio  tlao  is 


groator  than  tho  Ion  tlao  aoalo  (tg  »  upJ).  On  tho  othor  hand,  tho  tlao 
ooolo  for  tho  rolotlviatio  aaoa  lnoroooo  lo  tR  -  e*1  whoro  m  lo  tho  radiation 
frequenoy  and  la  aaouaod  ouch  groator  than  tho  oloctron  plaoao  froquoncy. 

Previous  roaoaroh  on  self-focusing  woo  prlnarily  ooncomod  with  tho 
pondoroaotlvo  self-focusing  effect.1®”1^  Max1 ^  presented  steady-state 
solutions  for  tho  radiation  bean  structure  due  to  ponderoootive  self-focusing 
In  a  nonrelativlstle  plaaaa.  In  this  treataent  it  was  assuaod  that  the  plasaa 
had  attained  equilibrium  with  the  radiation  field  and,  hence,  the  plasma 
density  was  represented  by  a  Boltzman-type  response,  n  -  exp  (ea  /T  ),  where 

P  9 

a  Is  the  ponderomotlve  potential  and  T_  is  the  electron  temperature. 

P  9 

Felber  ^  performed  a  similar  treatment  of  ponderomotlve  self -focusing  with  the 
lnolusion  of  relativistic  effects  and,  again,  the  assumption  of  an  equilibrium 


’f  T* 
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MMU;  distribution  mm  used.  Suoh  tftittwti,  using  sn 


g^lty 


•  srs  invalid  far  short  rsdlation  pulsss. 
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«  <  WGy.  Oh  sash  Mrt  tlae  stslsa.  ttea  platan  ions  oan  not  rs spond  to  ths 
fMwllft  far  on. 

fin  protlea  of  rolativlstlo  self -focusing  along,  valid  for  short  tiass 
before  Mas  ion  density  rasponds  to  pondsroaotivs  foros,  has  boon  oonsidsrsd  by 
Uni  at  al.,1*  Sprangls  and  Taag1*  and  Sohaldt  and  Horton.20  Muaarioal 


aiaalatlona  hava  alto  boon  attsaptsd  for  arbitrary  tints  including  ths  ooablnsd 
sffsots  of  relativist io,  pondsroaotivs  and  thsraal  sslf -focusing. 21 ,22 


Thin  problem  of  sslf-foousing  in  ths  transition  period  during  whloh  ths  plasna 
profile  begina  to  respond  to  pondsroaotivs  forest  and  thsraal  hsatlng  will  ba 
addrasssd  in  a  later  paper. 

Ths  analysis  presented  hers  will  be  oonoerned  only  with  relativistic 
sslf-foousing  on  a  tins  scale  sufficiently  short  so  that  ths  plasaa  density 
profile  dose  not  evolve  signlfloantly  under  ths  influence  of  the  radiation 
beaa.  This  laplles  that  the  pulse  length  of  the  radiation  beaa,  t^,  oust  be 
short  ooapared  to  is  and,  of  course,  long  coaparsd  to  a  radiation  period,  tR. 

The  relatlvlstlo  self-focusing  sffsot  is  analyzed  for  a  helically 
polarized  beaa  propagating  in  the  z  direction.  The  beaa  is  assuaed  to  be 
axially  syaastrlo  with  respect  to  ths  z  axis  and  has  a  profile  which  Is  only  a 
funotion  of  r,z  and  t.  An  equation  which  describes  the  envelope  of  the 
radiation  beaa  Is  derived.  The  envelope  equation  includes  diffraction  effects 
as  well  as  relstlvlstlo  plasaa  effects  and  has  a  fora  similar  to  a  particle 
beaa  envelope  equation.2^* 2**  In  ths  present  radistion  beaa  envelope  equation, 
dlffraotlon  effects  are  aanlfestsd  through  a  term  which  is  equivalent  to  beam 
eel ttanoe  in  ths  particle  beaa  envelope  equation. 
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a  helically  polarlsad  radiation 


propagating  within  a  cold 


*  The  footer  potential  of  tho  field  la  taken  to  have  the 


A^r, s,t)  •  A<r, a)  (coa<kx-wt)eK-  ain  (kx-wt)sy), 


the  aoplltude  A(r,s)  la  a  slowly  varying  function  of  r  and  x,  and  the 


frequency  «  la 


to  be  ouch  greater  than  the  effective  plasm 


frequency.  The  approximate  local  dispersion  relation  associated  with  the 


field  la  (1)  la 


«  •  ok  ♦  (c2k2  ♦  «t2(r,x)  )/2ok, 


where  k^  *  (k*  ♦  kp  Is  the  transverse  wave  number,  k*  «  k*  and  w*(r,z) 
la  the  effeotlve  background  plasm  frequency.  For  purposes  of  this  discussion 
the  effeotlve  plasm  frequency  is  written  In  the  fora 


.  J5 — r, 

"o  1'jp(r’x> 


•pC.*) 


2  1/2 

where  •»  •  (tviel  n  /m  )  Is  the  ambient  plasm  frequency,  n(r,z)  is 

po  1  ’  o  o 

the  oodlfled  electron  density  of  the  plasm  due  to  the  excited  plasm  wave, 
Y^(r,x)  •  (l*a2(r,z))1/2  is  the  relativistic  mss  factor  and  a(r,z)  - 
|e|A(r,x)/m0c2  is  the  normlized  laser  field  amplitude.  The  relativistic 
factor  arises  from  the  plasm  electrons'  relativistic  mss  change  due  to 
their  transverse  oscillations  in  the  radiation  field. 


>>«*■**«"  H«  tnyXop#  twtlon 

\i"  im  tU«  iMUon  an  envelope  aquation  la  derived  which  describes  the 

rtyn— I  in  of  the  radiation  beam  envelope  as  It  propagates  through  a 
The  derivation  la  aiailar  to  that  uaed  to  obtain  a  particle  beam 
envelope  equation211  in  the  aenae  that  the  ray  equations  in  (5a,  b)  have  a 
fora  which  la  aiailar  to  partiole  orbits  in  the  paraxial  approximation. 

To  this  end,  various  moments  are  taken  of  the  ray  equations  given  by 
(5).  Multiplying  (5a)  by  x  and  dx/dt,  (5b)  by  y  and  dy/dt  and  combining, 
yields  the  following  virial  and  energy  equations 


5  -  v2  ♦  fl2(r,z)r2  -  0, 

2  dtz 


(6a) 


g  *  -2<^>  #  ■  «. 

dt* 


(6b) 


where  r2  •  x2  ♦  y2  and  v2  •  (dx/dt)2  ♦  (dy/dt)2.  Substituting  (6a)  into  (6b) 
yields 


liE.fcttV,  *o2#"0. 

2  dt3  dt  dt 


(7) 


To  obtain  the  envelope  equation,  consider  a  thin  slice  of  the  beam  lying  in 
the  transverse.  The  transverse  segment  of  the  beam  is  assumed  to  consist  of  N 
individual  rays  traveling  through  it.  The  mean  square  radius  of  the  radiation 
beam  envelope  is  then  defined  as 


N 

R2(z, t)  -  <r2>  -  ij  l  rf 
M  i-1  1 


(8) 


6 


ti mr*  Pj  is  the  radial  position  of  the  ith  ray  and  <  >  donotas  tha  averaging 

operation.  Mathematically,  this  averaging  can  ba  written  in  taraa  of  an 
integral  over  tha  distribution  function  of  initial  ray  positions  at  tha 
initial  time  t0>f(?0,t0). 

<Q(r)>  -  /  dr  r  f(r  ,t  )Q(r) ,  (9) 

•  OO  0  0 


whara  Q (r)  is  a  quantity  dapandant  on  tha  ray  position  r  -  r(rQ,t). 
Performing  this  average  over  all  tha  rays  on  (7)  and  using  (8)  gives 


1 

2 


3t  <fl2f2>  * 


2  dr2 

<fl  3t> 


0. 


(10) 


In  order  to  express  tha  averages  in  (tO)  in  a  more  manageable  fora,  the 
individual  radial  ray  velocity  is  written  as  a  sun  of  a  mean  velocity  term  and 
a  residual  velocity  term.  That  Is  tha  individual  ray  velocity  is  written  as 


01) 


where  (r/R)dR/dt  is  defined  as  the  mean  radial  velocity  at  the  ray  position  r 
and  6vr  is  the  residual  radial  velocity.  Substituting  (11)  into  the  last  term 
in  (10)  yields 


<Q2  -  2<Q2r2> 

Qt 


.  2<o2?j;  >. 

H  P 


(12) 


To  further  simplify  (12),  notice  that 


•  <TT>-  I  ^  ♦  <fiT> 


03) 


l«Mt|  <r*v>  -  0.  Similarly,  assuming  <rn>  -  ©nIn  where  cn  it  a  constant. 


1  d_ 
n  dt 


<r*> 


,-n“1  *  „  .n-1 

<p  »>  •  c_R 
n 


♦  <rn-15v>. 

Qv 


(U) 


Bence,  the  conclusion  is  roach ad  that  <r*6v>  -  0  for  any  positive  integer  a. 

2 

Assuming  that  Q  (r)  can  bs  expanded  in  a  power  series  of  r  implies  that 
<02rlr>  -  0.  Thus,  (12)  becomes 


<02  £i> 

“  dt 


2<q2?2>  1  ** 
r  ?  R  df 


(15) 


Substituting  (15)  into  (10)  yields 


5  ^"T"  +  <q2''2>  +  2<Q2r2> 

2  dt3  dt  R 


(16) 


•2  2  2^2 

Since  the  last  two  taros  in  (16)  combine  to  give  R  d(R  <Q  r  >)/dt,  (16)  can 
be  written  as 


!_  rR2  dV 
dtl  dt2 


1  rdR2 

2  ldt 


2 


0. 


(17) 


•  •  2  2 
Integrating  (17)  and  taking  the  constant  of  integration  to  be  2c  e  gives 


0 


(18) 


d2!  ^  <K^r2>  _  cf 
d?  *  R3 

j  2  2 

uhsrs  •  S  /o  and  z  •  ct. 

Equation  (18)  describes  the  evolution  of  the  radiation  beaa  envelope,  R , 

ae  function  of  z.  Diffraction  effects  which  result  in  a  spreading  of  the  beam 

are  oontained  in  the  constant  of  integration  term  c.  The  form  of  the  envelope 

equation  for  the  radiation  bean,  (18),  is  similar  to  the  usual  particle  beam 

envelope  equation.  The  second  term  in  (18)  can  be  either  focusing  or 

2 

defoouaing  depending  on  the  sign  of  Q  .  In  the  present  situation,  because  of 

the  relativistic  transverse  oscillations  of  the  background  plasma  electrons, 

this  tern  results  in  focusing  of  the  radiation  beam.  The  third  term  in  (18) 

is  always  a  defoouaing  term  and  is  due  to  diffraction  effects. 

2 

The  diffraction  constant  e  can  be  evaluated  for  the  special  case  of  a 
Gaussian  ray  distribution.  Using  the  definition  in  (9)  and  requiring  <1>  -  1 
and  <r2>  •  R2  then  implies 


2  ~^q2/r2 

f(ro’V  -  j?  *  ' 


(19) 


Since  the  density  of  rays  is  proportional  to  the  intensity  of  the  radiation 
field,  this  implies  that  f(ro,tQ)  is  proportional  to  a  (rQ,t).  In  standard 
notation,  the  vector  potential  of  a  Gaussian  radiation  field  is  written  as 

2  2 

a  R  -r /R 

a(r*z) "  riff e  s’  (20) 

9 


where  R9(z)  is  the  spot  size,  Rgo  the  minimum  spot  size  and  aQ  is  the 
normalized  peak  radiation  amplitude.  In  a  vacuum,  Rs(z)  is  given  by 


9 


where  »  kR^/2  la  tha  Rayleigh  length25  and  z0  la  the  axial  location  of  the 

■tftijw  spot  size. 

•  2  • 

Since  f(r0, tQ)  la  proportional  to  a  (rQ,t),  then  by  comparing  (19)  and 

(20)  implies  that  the  mean-square  radial  ray  position  is  related  to  the  spot 

size  by  R  -  R  /2.  Rewriting  the  envelope  equation  (18)  in  terms  of  the  spot 
8 

size  Rs  gives 


2<?2K2> 


0. 


(22) 


The  integration  constant  eQ  can  now  be  evaluated  by  using  the  vacuum  solution 

(21).  In  vacuum,  <r2K2>  «  0  and,  hence,  c2  -  4/k2.  Notice  that 

-  A/*  -  R _ 0 .,  where  A  -  2*/k  and  8.  -  A/(trR _ )  is  the  well-known 

o  so  d  d  so 

diffraction  angle.25 


Solution  of  Envelope  Equation 

In  order  to  solve  the  envelope  equation  (22)  it  is  necessary  first  to 

-2  2 

evaluate  the  average  of  the  quantity  r  K  ,  where 


(1 


n(r,z)/n 
♦  a2(r,z)) 


1/2  ’ 


(23) 


where  f-?(?0,t)  and  z  *  ct  ♦  zQ. 

Reoall  that  the  averaging  can  be  written  in  terms  of  the  distribution  of 
initial  ray  positions  f(ro,tQ)  aa  ln  (9).  Notice  that  it  is  possible  to 
define  a  Eulerian  ray  .distribution  function  at  position  r  and  time  t  as 


follows: 


•1“  '-■••’  V-:;. 


.._,.  ,-^'.f  r:-  V; >,;*•. 


f(r, t)  -  /  dr0f(r0,t)«(r  -  r(?o,t)J. 


Consider  now  the  average  of  the  quantity  Q(r)  defined  by 


<Q(r)>  •  /  drrf(r,t)Q(r) 

o 


-  /  drr  J  dr0f(?0,t)«(r  -  r(ra,t))Q(r) 


o  o 


/  f<r  .t)Q(r) 
'  0  0  o 


<Q(r)>. 


Hence,  the  averaging  <Q(r)>  can  be  replaced  by  an  integration  over  the 
Euler lan  spatial  coordinate  r  providing  the  distribution  f(r,t)  is  known. 

Since  there  la  no  hope  in  determining  an  analytic  expression  for  the 
distribution  f(r,t),  in  order  to  obtain  an  estlaate  of  <r2K2>,  it  will  be 
asauaed  that  the  radiation  beaa  retains  a  Gaussian  profile  as  it  propagates 
through  the  plasaa.  Thus,  f(r,t)  will  be  approxiaated  as 


f(r,z)  -  -r-.-  a2(r,z). 


where  z  •  ot  ♦  z0  and  a(r,z)  Is  given  by  (20).  Notice  that  <1>  -  1  and 
<r*>  •  R2/2.  Using  this  expression,  then  <r^K^>  becomes 

9  9  2  •  _  -  «  n(r,z)/n, 

<r  A  *2(^t)  l  * r  •  <r-i)  fe 


2  •  __  n(s)/nrt 

’^ck2^  1  ds(1-s)e  - •2~•3~T/2,  ' 

c*  o  (Hx  ze  V  2 


11 


~  9  9  90  O’ 

9&P  the  ossa  of  9  constant  density  background  plasma,  n(r)  -  nQ  for 
'¥  4  '-HumM  Mr)  •  0  for  r  >  Rp,  then  (27)  oan  be  evaluated  analytically.  For 
this  ease,  letting  y  -  e“*,  then 


te2 

-  -§jj|  x2  {h(x)  -  g(x)(1-lnyp)  -  }  •  (28) 

Nfcere  h  -  [( 1*x"2)1/2-1 ]/2,  g  -  [0+x“2y  )1/2-l  ]/2  and  yn  -  exp(-2R2/R2). 

P  P  p  O 

It  la  insightful  to  write  the  envelope  equation  in  terms  of  an  effective 
particle  located  at  R#(z)  moving  in  a  potential  V(Ra).  Equation  (22)  becomes 


ttiere 


av(R#) 

ML 


9 


[h(x)  -  g(x)(1-lny&)  -  -  -y-^, 

k  R_ 


(29) 


(30) 


with  x  -V(aoRso)  and  yp  *  «*P(-2R2/R2). 


Envelope  Behavior 

To  analyze  the  behavior  of  the  radiation  envelope,  it  is  most  convenient 

to  use  the  normalized  envelope  radius  x  ■  R  /(a  R _ ) .  In  the  limit  R  ♦  •, 

s  o  so  p 

then  (30)  can  be  written  as 


(3D 


where 


f  *  x*3  -  16ox[((1  ♦  x~2)1/2-l)  ♦  2tn2 


-  21n((1*x‘2)1/2  *  1 )],  (32) 


»o  -  (»>/<kJlJ.*»2  «■*  <■  -  (v*0R«/(*c»2- 


The  above  aquation  describes  the  position  of  a  particle  x(t)  moving  in  the 
effeotive  potential  V(x,a).  In  the  expression  for  3V/3x,  the  first  term  on 
the  right  of  (32)  represents  vacuum  diffraction  whereas  the  term  proportional 
to  a  represents  the  relativistic  self-focusing  of  the  plasma. 

It  is  interesting  to  note  that  the  shape  of  the  potential  V(x)  depends 
only  on  the  parameter  a.  As  will  be  shown  below,  a  can  be  written  as  the 
laser  power  over  the  critical  power.  For  a  >  1,  the  potential  V(x)  has  a 
minimum  and  bounded  oscillatory  solutions  for  x(t)  are  possible. 

Expanding  V(x)  for  large  x  (small  aQ)  gives 


H  -  »-3  -  ox"3 


3x 


(33) 


For  o»1,  there  are  no  net  diffraction  or  focusing  forces  (to  leading  order) 

and  a  matched  radiation  beam  is  possible.  In  this  limit  (x  »  1),  the 

envelope  diffracts  for  a  <  1  and  focuses  for  a  >  1 .  Since  the  total  radiation 

2  2  2 

power  is  given  by  P  ■  (m  e  «a  R  )  /(8ce  ),  the  parameter  a  can  be  written 

6  O  SO 

as  o  -  where  the  critical  power  is  given  by 


crit 


no2  2  2 

*=(-!-)  ^ 

u 

po 


17.4  x  10 


9  « 


Watts. 


(34) 
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•  >v  I  ’  ,<§L"’'CSj‘ 


*,  •;«  ■■ 

;  ‘f  '  •  •  • 


I'Ssl*#*’"*  *  «  *  '«'*»• 


:  -  «ph  w  *v  •• 

j  .,  /=* 

til)  indicates  that  for  x  »  1 ,  and  a  >  1  *  than  the  radiation 


« 

■a 


ffcMtti  until  x  »  1  and  (33)  is  no  longer  valid.  Expanding  (32)  for  small 

Sfe  s 

Jgfaf ,«;j  gWb 


m 


m . 

«*C?  -  '■ 


*  x~^  -  16a. 


Equation  (35)  atatea  that  for  a  fixed  a  >  1 »  focusing  will  continue  until  the 
dlSfraotion  term  dominates  when  x^  <  1/(16a).  Hence,  at  a  sufficiently 
alMll  x,  the  envelope  will  be  reflected  back  out  towards  its  original  width. 
"'Ulan  this  occurs  the  envelope  radius  x  will  either  oscillate  between  its 
minimum  value  at  reflection  and  its  original  value,  or  it  will  continue  to 
diffract  indefinitely,  depending  on  the  envelope's  initial  "velocity"  (the 
initial  slope  in  x(t)).  For  a  given  a  >  1  and  an  initial  x  »  1,  if  the 
initial  slope  dx/dt  is  sufficiently  small,  then  x(t)  will  oscillate  between 
its  initial  value  and  its  value  at  reflection.  If  dx/dt  is  initially  large, 
x(t)  will  initially  decrease  to  its  minimum  value  at  reflection  and  then 
Increase  indefinitely. 

The  exact  shape  of  V(x)  for  a  given  a  must  be  determined  numerically.  In 
general,  if  a  >  1  then  there  will  exist  a  finite  well  whose  minimum  occurs 
at  Xf.  As  a  increases,  the  depth  of  this  well  increases,  the  well  becomes 
narrower  and  the  location  of  the  minimum  xf  decreases.  The  potential  V(x) 
is  plotted  in  Fig.  1  for  a  •  5.0  and  Fig.  2  for  a  -  1.2.  Notice  that  xf  -  0.4 
for  a  -  5.0  and  x^  -  1.8  for  a  *  1.2. 

Recall  from  the  vacuum  solution  that  the  minimum  of  Rs/Rgo  in  vacuum  is 
unity.  Hence,  the  minimum  for  x  in  vacuum  is  x  -  l/aQ.  Ideally,  it  is 
possible  to  have  a  nonoscillating  envelope  of  constant  radius  if  the  radiation 
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CM  pltMi  A  t  -  t0  iflth  x<t0)  -  xf  and  dx/dt  -  0.  This  corresponds 
to  the  radiation  anttrlng  tha  plaama  precisely  at  tha  bottom  of  tha  potantial 
mail  ulM  zero  velocity.  In  vacuum,  dx/dt  •  0  at  tha  minimum  spot  size 

»  Rm  and,  hence,  x  -  1/aQ.  In  order  to  achieve  a  nonoao mating  aolution, 
it  ie  required  that  x(t0>  •  xf  •  1 /aQ.  Kanoa,  for  a  -  5.0  (a  -  1.2)  a 
nonoao  mating  match  ad  beam  would  ba  poaalbla  for  a0  -  2.5  (a0  -  0.55) 
provided  the  radiation  anvalopa  an tar a  tha  plaama  with  zero  alope,  dR9/dz  -  0. 
Since  currant  technology  limlta  aQ  <  1,  operation  in  this  conatant  anvalopa 
node  la  poaalbla  only  for  a  alightly  greater  than  unity. 

In  practice,  it  may  not  ba  poaalbla  to  achieve  a  radiation  baaa  which 
entera  a  plaama  with  dRs/dz  ■  0  and  x  •  xf.  However,  it  may  ba  poaalbla  to 
operate  near  x  -  xf  with  a  small  initial  dx/dt.  In  thla  caaa,  whan 
x  •  x^  ♦  <x  with  6x  small,  tha  anvalopa  aquation  (3D  oan  ba  expanded  to  give 


d  .  2. 

“2  ,x  *  *  *o.4x- 
dt 


(36) 


where 


2  *¥o  ,  ,  ,3 

*0.  •  T  '  *“r  - 


(37) 

•f  (DxJ  ♦ 

For  tha  paramatara  a  ■  1.2,  aQ  -  0.55,  Rso  -  0.1  cm  and  w/unA  -  10,  this  gives 


po 


2  2  -6 

x„  •  1.8,  l  -  2wo/«.  -  43.5  cm  and  w„/.i  •  3  «  10  .  The  Rayleigh  length 

r  os  os  09  po 

for  thla  oasa  is  •  4.0  cm. 

Figures  3a>6  show  numerical  results  for  the  envelope  behavior  for  a  test 
oaaa  where  a  •  1.2,  Rso  -  0.1  cm  and  w/«^  -  10.0.  The  critical  power  is 


po 


.11 


given  by  Perlt  *  17  *  10  W,  the  well  minimum  occurs  at  xf  -  1.8  and  the 

oaolllating  wavelength  in  the  wall  is  11  44  cm.  Figure  3  shows  operation 

at  aQ  «  1/xf  *  0.55 *ahd  with  zero  initial  slope,  dRs/dz  -  0.  This  is  the 

1  2 

condition  for  a  nonoscillating  matched  beam.  Here  P  *  2.1  *  10  W, 


*  2.*  *  lO^sec”1,  m  •  2.4  *  lO^seo"*1  and  the  Rayleigh  langfch  la 
ig  *  4.0  oa.  Plgura  4  shows  a  oaaa  with  a0  slightly  lass  than  1/xf  *  0.55 
with  zaro  Initial  slope,  whila  Pig.  5  shows  a  casa  with  aQ  -  0.55  and  a  small 
initial  slopa.  In  both  thaaa  easas,  the  envelope  shows  a  snail  amplitude 
oscillation  whose  wavelength  is  given  by  x  •  1  *  44  on.  Figure  6  presents  a 

ease  where  a0  «  0.55  (initially,  x  »  1)  and  with  a  significant  initial 
slope.  Here  the  initial  slope  la  determined  by  specifying  zQ  where 
R 8(z)  -  R80(1  ♦  z‘/zR)  before  entering  the  plasaa  at  z  •  0.  For  this  case, 
the  radiation  initially  focuses  down  to  a  value  significantly  less  than  R90 
("over-focusing")  and  then  is  reflected  and  diffraots  indefinitely. 

When  larger  laser  powers  are  used  (larger  values  of  a)  the  results  are 
qualitatively  similar  except  the  potential  well  is  now  deeper,  narrower  and 
has  a  minimum  occurring  at  a  smaller  value  of  Xf.  Por  example,  when  a  -  5, 

Rso  “  *1  om  and  «/«po  "  10*  then  xf  *  0.4i  and  XQ#  •  8.9  cm.  A  nonoscillating 
matohed  beam  requires  ao  •  2.5  and  no  initial  slope.  For  a0  •  2.5,  this  gives 
*  1.7  *  1012W,  P  -  8.7  «  1012¥,  *  •  1.1  *  1 01 2sec"1 ,  «  -  1.1  «  lO^sec*1 

OTIC  pO 

and  zR  •  1.8  cm. 

numerical  Simulation 

In  order  to  oonflrm  the  above  results,  the  wave  equation  for  the 
radiation  is  solved  numerically.  For  a  plasma  with  constant  density,  and  a 
driving  current  that  results  from  the  relativistic  electron  oscillations,  the 
wave  equation  can  be  written  as 


Where  tM  slONty  vary  lag  amplitude  approximation  was  used.  This  aquation  la 
gal  rail  numerically  with  64  x  64  transverse  Fourlar  nodes. 

Tbs  laltlal  transverse  radiation  pronia  was  a  Gaussian  with  zero 
wavefront  ourvature.  As  the  radiation  propagated  through  the  plasna,  the 
profile  dors loped  finite  but  saall  aaplltude  "skirts"  indioatlng  higher  order 
transverse  Gaussian  nodes  are  involved.  The  1/e  radius  (at  which  the 
amplitude  is  1/e  that  on  axis)  is  chosen  to  be  the  aeasure  of  the  radiation 
envelope.  This  is  plotted  as  a  function  of  z  in  Fig.  7  for  the  case  when  the 
initial  nornallzed  envelope  x  ■Ra/(Rsoa0)  is  somewhat  greater  than  that  needed 
for  a  matched  bean.  Figure  7  indicates  that  the  envelope  oscillates  about  the 
predicted  matched  beam  radius.  The  oscillatory  behavior  is  no  longer  simply 
harmonic,  as  implied  by  the  above  theory,  and  the  oscillation  wavelength  is 
approximately  twice  the  predicted  value.  Notice  that  after  the  initial 
transient,  the  maximum  of  the  envelope  oscillation  is  only  about  2/3  of  the 
initial  radius  of  the  radiation.  This  discrepancy  is  probably  due  to  the  fact 
that  the  profile  deviates  from  its  initially  Gaussian  shape.  For  the  case 
where  the  value  of  the  initial  normalized  envelope  x  approaches  that  required 
for  a  matohed  beam,  the  oscillation  amplitude  becomes  very  small  and  the 
oscillatory  behavior  becomes  more  erratic. 

In  general,  the  numerical  simulations  support  the  conclusion  that  the 
radiation  envelope  oscillates  about  the  value  required  for  matched  beam 
propagation. 
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Tfet  <Wti  results  indioate  that  relativistic  electron  quiver  motion  leads 
to  eodsaoed  focusing  of  radiation  beams.  Physically,  relativistic  quivering 
ef  Uie  electrons  leads  to  an  effective  decrease  in  the  plasma  density  via  the 
electron  pis sms  frequency,  -  nfr.zjw^/tn^)  where  Y  •  1  ♦  a  (rtx).  A 

radiation  field  a2(r,z)  peaked  on  axis  at  r  -  0  will  produce  an  index  of 
re fr motion  profile  n(r,z)  peaked  on  axis  such  that  3n/3r  <  0,  where 
q  -  ck/e  •  (1  -  w‘(r,z)/w  )  .  Since  3q/3r  <  0,  this  implies  focusing  of  the 

radiation  beam. 

The  effects  of  relativistic  self-focusing  on  the  development  of  the 
radiation  envelope  Rs(z)  are  best  understood  through  the  analogy  or  a  single 
particle  with  orbit  R#(z)  moving  in  an  effective  potential  V(Ra).  As 
dlsoussed  above,  the  shape  of  the  potential  V  is  determined  by  a  single 

a 

parameter  a  -  p/p0Plt»  where  the  critical  power  is  given  by  pcrlt  •  17  *  10 

(m/*P©)2tfatts.  Provided  a  >  1  (laser  power  greater  than  the  critical  power), 

then  there  exists  a  minimum  in  the  effective  potential  V  located  at 

x.  •  xAa),  where  x  •  R  /(R  a_).  As  a  is  increased,  the  well  depth 
r  r  «  so  o 

increases  and  becomes  narrower,  and  the  location  of  the  minimum  decreases. 
The  existence  of  such  a  well  implies  that  bounded  solutions  for  Rs(z)  are 
possible  where  Ra(z)  oscillates  between  the  two  reflection  points  associated 
with  the  effective  potential.  For  example,  if  a  >  1,  a  radiation  beam 
entering  the  plasma  with  an  initial  value  of  x  greater  than  xf  will  initially 
begin  to  focus  (x  decreases).  When  this  occurs,  two  outcomes  are  possible: 

1)  the  normalized  envelope  x(z)  will  continue  to  decrease  to  some  minimum 
value  at  which  it  will  be  reflected  and  expand  indefinitely,  or  ii)  the 
envelope  x(z)  will  remain  bounded,  oscillating  Indefinitely  between  the 
minimum  reflection  point  and  its  initial  value.  The  occurrence  of  one  or  the 
other  of  the  above  two  cases  depends  on  the  initial  "convergence  angle"  of  the 


envelop*  x(z)  «s  it  enters  the  plasma  or,  more  precisely,  on  the  Initial  slope 
dltg/dz.  If  dRa/dz  is  too  large  then  case  i)  occurs:  the  envelope  is 
reflected  and  expands  indefinitely.  Zf  the  initial  dRa/dz  is  sufficiently 
small,  then  case  11)  occurs:  the  envelope  becomes  trapped  in  the  potential 
well  and  oaoillates  about  the  minimum  xf.  Under  a  special  set  of  initial 
oonditlons,  dRa/dz  •  0  and  x  -  x^,  it  is  possible  to  have  a  matched  beam  with 
a  constant  envelope. 

The  following  approximations  were  made  in  this  analysis:  the 

mathematical  model  used  in  the  above  calculations  was  based  on  the  ray 

equations  of  geometrical  optics.  These  ray  equations  then  used  a  dispersion 

relation  for  an  electromagnetic  wave  in  an  unmagnetized  plasma  in  which  it  was 
2  2  2  2 

assumed  k  »  and  a  »  The  radiation  beam  was  initially  assumed  to  be 

Gaussian  and,  furthermore,  it  was  assumed  to  remain  Gaussian  as  it  propagated 

through  the  plasma,  a(r,z)  -  (a0Rao/Ra)exp(-r2/R^(z)).  The  last 

approximation,  which  restricts  the  analysis  to  short  pulses,  assumes  that  the 

background  plasma  density  remains  constant  and  does  not  evolve  under  the 

influence  of  the  transverse  ponderomotlve  force.  Such  an  approximation  should 

be  valid  for  short  pulse  times  t  <  R  /C  . 

o  s 

The  results  discussed  above  are  in  qualitative  disagreement  with  those  of 
Schmidt  and  Horton,20  who  also  analyzed  relativistic  self-focusing  for  short 
pulses.  They  found  a  similar  expression  for  Pcplt,  beyond  which  self-focusing 
occurs.  They  claimed,  however,  that  for  P  >  Pcrit  the  radiation  beam  would 
collapse  down  to  some  minimum  radius  and  neither  oscillate  nor  reflect.  This 
is  similar  to  the  results  of  Max  et  al.,18  in  which  it  was  determined  that 
self-focusing  occurs  when  a  certain  threshold  is  surpassed.  Their  analysis, 
however,  did  not  predict  the  envelope  behavior  once  self-focusing  occurs.  The 
above  numerical  simulation  of  the  wave  equation  Indicates  that  the  correct 


envelope  behavior  is  an  oscillation  about  the  value  required  for  a  matched 

beast. 

The  results  of  Felber,1®  which  included  relativlstio  focusing  along  with 
that  of  the  equilibrium  ponderoaotive  effect,  are  in  qualitative  agreement 
with  those  presented  above.  Specifically,  Felber16  found  that  the  behavior  of 
the  radiation  envelope  could  also  be  described  as  a  particle  in  an  effective 
potential.  Provided  the  laser  power  was  sufficiently  high,  the  potential 
exhibited  a  minimum  and,  hence,  either  osoillating,  bounded  solutions  or 
reflected,  diverging  solutions  were  possible.  The  critical  power  in  Felber's 
case  was  similar  to  that  found  by  Max,15  pCrit  *  2  *  lO^KevJw2/^  Watts, 
which  is  typically  less  than  that  disoussed  above  when  the  effects  of 
ponderomotlve  focusing  are  neglected.  The  main  limitation  with  the  analysis 
of  Felber*®  and  Max1^  is  that  both  used  an  equilibrium  density  response, 
n  -  exp(e+p/T)  where  is  the  ponderomotlve  potential.  Such  an  equilibrium 
density  response  is  invalid  for  the  leading  segment  of  the  radiation  pulse. 
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Figure  2.  The  effective  potential  V  (x)  as  a  function  of  x 
for  o  •  1.2.  The  well  minimum  occurs  at  x-  -  1.8 


Figure  3.  The  radiation  envelop*  Ra(z)  for  the  paraaeters  a  -  1.2 
II  -  0.1  ca,  m/u _ ■  10.0  and  a„  -  0.55.  For  this  run 

•°  it  po  ?3  -1 

pcrit  *  H  and  e  -  2.*x10  ■’sec  .  Initially, 

x  •  x„  ■  1 .8  and  dx/dz  -  0. 


0 

Figure  4 
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The  radiation  envelope  Ra(z)  for  the  parameters  a  -  1.2, 

RgQ  -  0.1  cm,  «/«p0  ■  10.0  and  aQ  -  0.54.  For  this  run 

^orit  *  ITxIO^W,  a  -  2.4x1 01  ^sec*1 ,  zR  •  4.1  cm  and 

I..  *  44  cm.  Initially,  x  »  1.02  x.  and  dx/dz  -  0. 
oo  r 
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Figure  7.  The  1/e  radial  width  In  terna  of  x  as  a  function  of 

z  (In  oa)  aa  determined  by  numerical  aimulation  of  Eq.  (38).  Fi 
thia  run,  o  •  5,  u/id _ -  10,  a  * 

nn  n 


2.21  and  R 


0.11  cm. 
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